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Abstract 

The Kamioka Liquid scintillator Anti-Neutrino Detector (KamLAND) is sensitive to the neutrino 
event spectrum from (mainly Japanese) nuclear reactors in both the energy domain and the time 
domain. While the energy spectrum of KamLAND events allows the determination of the neutrino 
oscillation parameters, the time spectrum can be used to monitor known and unknown neutrino 
sources. By using available monthly-binned data on event-by-event energies in KamLAND and 
on reactor powers in Japan, we perform a likelihood analysis of the neutrino event spectra in 
energy and time, and find significant indications in favor of time variations of the known reactor 
sources, as compared with the hypothetical case of constant reactor neutrino flux. We also find that 
the KamLAND data place interesting upper limits on the power of a speculative nuclear reactor 
operating in the Earth's core (the so-called georeactor); such limits are strengthened by including 
solar neutrino constraints on the neutrino mass and mixing parameters. Our results corroborate 
the standard interpretation of the KamLAND signal as due to oscillating neutrinos from known 
reactor sources. 

PACS numbers: 14.60.Pq, 28.50.Hw, 26.65. -ft, 91. 35.-x 
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I. INTRODUCTION 



The Kamioka Liquid scintillator Anti-Neutrino Detector (KamLAND) 0, 0| is sensitive 
to oscillations 0,111 of reactor neutrinos [5] over longbaselines ((L) ~ 180 km). The neutrino 
disappearance effect observed in KamLAND la [7|, |8j provides an independent confirmation 
of the matter-enhanced adiabatic solution jij [h| to the solar neutrino problem 11 , [T~ 



Il8l | at large mixing an gle (LMA), with best-fit oscillation parameters 
(5m 2 , sin 2 12 ) ~ (7.9 x 10" 5 eV 2 ,0.31) 0,13 in standard notation In addition, the 

current KamLAND statistics and energy resolution allow to track the oscillatory pattern of 
reactor neutrinos in the energy domain for about half a period • 

Being a real-time detector, KamLAND can also track neutrino source variations in the 
time domain. In particular, significant power variations of some Japanese reactors occurred 
during data taking 0] , leading to expected variations in the KamLAND neutrino event rate 
0. The KamLAND sensitivity to time variations was estimated to reach potentially the 
~ 2.3<T level through the unbinned test proposed in where only time information (and 
no energy information) was considered. 

There is also, in principle, an interesting interplay between time and energy information 
in KamLAND. In the presence of neutrino oscillations, time variations of reactors placed 
at different distances produce time variations of the energy spectrum. Figure 1 shows, e.g., 
that by "switching off" one of the most powerful nuclear reactor plants in Japan (namely, 
Kashiwazaki) one gets not only an overall decrease of the spectrum normalization, but also 
a slight displacement of the dip in the oscillatory pattern. In the same figure, one can also 
see the effect of a hypothetical reactor at the center of the Earth (the so-called georeactor 
[21]), which would increase the KamLAND spectrum by a factor which is constant in time 
but, in general, not uniform in energy. A joint analysis in the energy and time domain would 
be appropriate to study such effects. 

So far, the KamLAND collaboration has published only one test of the time-variation 
hypothesis, which makes use of a relatively coarse time binning and of no energy information. 
The results are shown in the first figure of Ref . , where the observed event rates — grouped 
in five data points — are plotted against the unoscillated reactor neutrino flux, and a positive 
(expected) correlation is seen to emerge. However, the statistical difference between the two 
extreme cases in this test (with and without time variations of the neutrino flux) is only 
Ax 2 = 3.3 0, i.e., smaller than 2a. At a similar significance level, the extrapolation to 
"zero reactor power" is consistent with the known background, but yields poor constraints 
on possible unidentified sources such as the georeactor 

The power of a time-variation test — as the one described in — can be improved by 
exploiting additional information. For instance, daily data about individual Japanese reactor 
operations are available to the KamLAND collaboration, through an agreement with the 
power companies ji^. In principle, these data allow one to perform detailed likelihood 



analyses of the event spectra not only in the energy domain (as those, e.g., in 0, 0, EH 
IIll^Ql2ll2El^l2ll^Sl2ill2ll23)but also in the time domain, thus providing 
statistically more powerful tests of reactor power variations and of the georeactor hypothesis. 
Unfortunately, daily reactor data are classified [22] |. 

Recently, monthly-binned data from nuclear reactors and from KamLAND have become 
publicly available. In particular, average Japanese reactor powers in each calendar month 



can be found at |3J|. The sequence of published KamLAND events |8( in monthly bins, 



tmy 

together with the corresponding detector livetime (in seconds), can be found in [35(. The 
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availability of these data has prompted us to extend the likelihood analysis of KamLAND 
data in the energy domain (event by event) [f| 0, |25[ so as to include the time domain 
(monthly binned). We find that the joint maximum-likelihood analysis in energy and time 
can provide a significant (~ 3cr) indication in favor of time variations of reactor powers, as 
compared with the case of average constant powers. In addition, we find no indication in 
favor of a georeactor contribution, and we set upper bounds on its power. In both cases, 
we discuss the role of additional solar neutrino constraints on (5m 2 , sin 2 6 I 12 ). Our results 
corroborate the standard interpretation of the KamLAND signal as due to flavor oscillations 
of neutrinos coming from known reactor sources. 

The structure of this paper is as follows. In Sec. II we reproduce, as a preliminary but 
relevant check, the official KamLAND unbinned likelihood analysis in the energy domain 
In sec. Ill we extend the analysis to the (monthly binned) time domain, and show that 
significant indications in favor of reactor time variations emerge from the data. In Sec. IV 
we discuss the effects of a hypothetical georeactor, and set upper bounds on its power. We 
summarize our results in Sec. V. 

A final remark is in order. Our results, although encouraging, cannot — and must not — be 
taken as a substitute for future, official KamLAND tests of hypotheses about the reactor 
sources. In fact, as described in the following, our approach involves some unavoidable 
approximations, which could be easily removed by the KamLAND collaboration — possibly 
leading to somewhat different results. Nevertheless, we think that our approximate analysis 
in the energy-time domain may represent an interesting step beyond previous KamLAND 
data analyses, where the time information is absent. 



II. LIKELIHOOD ANALYSIS IN ENERGY 

The KamLAND experiment has collected so far A Q b s = 258 events in a fiducial mass 
M = 543.7 tons, during a total livetime At = 515.1 days Q: 

M ■ At = 0.766 kTy. (1) 

Details on the likelihood analysis of the energy spectrum of such events are available at 
0]. In this Section we reproduce the results of the official KamLAND likelihood analysis in 
energy 0, before generalizing it to the time domain in Sec. III. 

In general, the KamLAND unbinned likelihood function C can be written as 0, 0, : 



£ ^rate X -£ s hape X /3 S y S t , (2) 

where the three factors embed information on the total event rate, on the spectrum shape, 
and on systematic uncertainties. The evaluation of £ implies a detailed calculation of the 
absolute spectrum of events (signal plus background), whose ingredients are briefly described 
below. 



A. Reactor input 



The reactor signal in KamLAND is essentially generated by 20 nuclear reactor power 
plants (16 in Japan and 4 in Korea) located at different distances Lj and characterized by 
different thermal powers Pj h [2]. For Japanese reactors (j = 1, . . . , 16), the sequence of 



3 



monthly-averaged thermal powers Pj^ (where m is a monthly index) can be recovered from 
the corresponding sequence of average electric powers Pj^ available in j34|, by using the 
relation P th ~ 3 P el |36| . The time interval of interest for the current KamLAND analysis 
spans m = 1, . . . , 23 months, from March 2002 to January 2004 included 0, H^]. In each 
month, the KamLAND detector livetime At m (with J2m — At) is given in 35|| . The 



average thermal power of the j-th Japanese reactor during the total KamLAND livetime At 
can thus be approximately estimated as 



y>26 pth Af 

pth ^ l^m=l r jm L ^ i m , - 

^ " At ' {6) 



where we are implicitly neglecting variations of the reactor powers (and of the detector 
livetime) over time scales shorter than a month. 1 We have not found monthly information 
about the four Korean reactor plants (j = 17, ...,20), which we simply assume to have 
constant powers (Pj^ = Pf 1 ); where Pj h is taken as a typical fraction (80%) of the nominal 
thermal power quoted in |2j. 2 For all reactors, the average fuel components qj (/ = 1, . . . , 4) 



are taken as 

235tj . 238TJ . 239p u . 241p u = q ggg . q q7q . q^qj . Q m? ^ 

at all times, with average fission energies Ef = 201.7, 210.0, 205.0, and 212.4 MeV, respec- 
tively [3?}. We do not have enough information to implement fuel burn-up corrections 
to individual reactors. 

Within the above approximations, the time-averaged differential neutrino flux at Kam- 
LAND (number of neutrinos per unit of time, area, and energy) is then given by [H[ 

dE v fy^AirtfEfdE,, ' 1 ' 

where we assume, for the /-th spectral component, the parametrization |3 

^ = exp(^ + aK + a^ 2 ) , (6) 

the a{ coefficients being reported in [3^]. In the presence of oscillations, each j-th reac- 
tor term in Eq. (jHJ) must be multiplied by the corresponding neutrino survival probability 

Pee{E u , Lj). 

We have made two reassuring checks of the above reactor power input. As a first check, 
we have estimated the total integrated thermal power flux over the detector livetime, 

£ 3 ^ = 697 J/cm\ (7) 
in good agreeement with the official KamLAND value of 701 J/cm 2 Q. 



1 This approximation could be removed by using, e.g., daily data, which are available only within the 
KamLAND collaboration. 

2 This is a minor approximation, since Korean reactors contribute only ~ 3% to the KamLAND signal. 
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As a second check, we have calculated the so-called integrated fission number flux 

and its (binned) distribution over the reactor distance L. Figure 2 shows that our results 
are very close to the official KamLAND ones, as taken from |8[. 



B. Detection input 

Given the differential neutrino flux in Eq. (0), the time-averaged energy spectrum of 
reactor events in KamLAND (number of expected events per unit of prompt positron energy 
E) is given by 

S(E) = enMAtjdE v ^-JdE> r{E, E') , (9) 

where e = 0.898 is the overall efficiency (after all cuts 0), n is the target density (0.848 x 10 29 
protons/ton) 0, HE,E') is the energy resolution function (with Gaussian width equal to 
7%(-E"/MeV) 1 / 2 ) Pi, and a is the inverse beta decay cross section, estimated as 

^%r^ ^ a{E v ) 8{E V — E' — 0.782 MeV) , (10) 

with a(E u ) taken from In r(E,E'), we allow for a systematic offset of the prompt 

(true) energy scale, 

E' — E'{1 + a) , (11) 

with standard deviation a a = 2 x 10~ 2 0. 

Above the current analysis threshold (E thr = 2.6 MeV), we estimate a total of 377.3 
reactor events in the absence of oscillations. This value is about 3% higher than the official 
KamLAND estimate (365.2 events 0); we obtain a +3% difference also in comparison with 
older data (89.7 events against the official 86.8 estimate 0). We have not been able to 
trace the source of this modest systematic difference, which we choose to compensate "ad 
hoc" in the following, through a fudge factor / = 0.97 multiplying the right hand side of 
Eq. ©. 3 

Finally, one must consider the background energy spectrum B(E) expected over the 
livetime At. This spectrum has three main components, as described in detail in 0,0: the 
accidental background B\, the 8 He- 9 Li background B2, and the 13 C(a,n) 16 background 
-B3. While Bi and B2 can be estimated with very small uncertainties [7| (that we set to 
zero), the normalization of the third background is poorly known in both its low-energy 
(< 5.4 MeV) and high-energy (> 5.4 MeV) components |g| (B' 3 and B 3 , respectively). We 
then assume free normalization factors (a' and a") for such components. In conclusion, we 
take the absolute background spectrum as 

B(E) = B^E) + B 2 (E) + a'B' 3 (E) + a"B' 3 '{E) , (12) 

where the B\, B2, B' 3 and B 3 components are taken from ,8|, while a' and a" are free 
(positive) parameters. 



3 This small adjustment (3%) is only ~ 1/2 of the KamLAND normalization error (6.5%). Removal of such 
adjustment does not appreciably change any of our results. 
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C. Likelihood function and oscillation parameters 



The absolute energy spectrum of events expected above the analysis threshold can always 
be factorized into the total number of events iVtheo times the probability distribution in 
energy D(E), namely 

S(E)+B(E) = N theo -D(E) (13) 

with 

f dED(E) = l. (14) 

We remind that both iVtheo and D(E) depend on the systematic energy offset a, as well 
as on the free background parameters a' and a". In the presence of oscillations, they also 
depend on the the mass-mixing parameters (5m 2 , sin 2 #12). 4 

Given the previous definitions, the first likelihood factor in Eq. (J2J) can be written as (see 
also |25|): 



rate 



2na 



exp 



rate 



1 / N theo (5m 2 , sin 2 9 12 ; a, a', a") - N ohs ^ 



"rate 



(15) 



where iV bs = 258 is the total number of observed events 0], and the total error is the sum 
of the statistical and systematic (s = 6.5% 0,I1|) uncertainties, 



a 



rate 



iVtheo + (S N< 



theo 1 



(16) 



The second likelihood factor in Eq. (J2J) is the product of the probability that the z-th event 
(i — 1, . . . , iV bs) occurs with the observed energy E i: 



258 



£shape = II D ( E i I ^ m2 ' Sin2 ^12 5 a : «'> a ") 



(17) 



where the energy set {E{\ is given in [8J]. The third and last likelihood factor in Eq. (J2J 
embeds the penalty for the systematic offset a in Eq. (jllj) . 



-■syst 



exp 



2 VOn 



(18) 



In general, further penalty factors could account for additional KamLAND systematics (not 
included here for lack of detailed published information). 
Finally, the standard x 2 function is obtained as 



X 2 (5m 2 , sin 2 6 12 ) = —2 In max C(5m 2 , sin 2 6*12 ; a, a', a") 

{a,a' ,a"} 



(19) 



Bounds on the oscillation parameters can be found by plotting isolines of the function 

A X 2 = X 2 ~ min X 2 ■ (20) 

{6m 2 ,sf 2 } 



4 In this work we do not consider subleading three- neutrino oscillation effects, i.e., we assume 613 = in 
standard notation. Within current bounds (sin 2 #13 < few%) we do not expect this app roximation to be 
crucial. We also neglect small Earth matter effects on reactor neutrino propagation |28j . 
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The values A X 2 = 4.61, 5.99, 9.21, and 11.83 correspond to 90, 95, 99, and 99.73% C.L. for 
two degrees of freedom. 

Figure 3 shows the bounds on the oscillation parameters from our likelihood analysis of 
the KamLAND energy spectrum. The confidence level isolines are in very good agreement 
with the official ones reported in Fig. 4(a) of 0, modulo the different scales chosen for the 
axes. 5 These results, together with the previous checks in this Section, demonstrate that we 
can reproduce, to a good accuracy, both the input and the output of the official KamLAND 
likelihood analysis of the energy spectrum. This check is also relevant to appreciate, in 
the next Section, the (small) differences induced by including the time information in the 
likelihood analysis. 



III. LIKELIHOOD ANALYSIS IN ENERGY AND TIME 



The information reported in [35[ allows to separate the global KamLAND set of 258 
event-by-event energies into 23 monthly subsets I m , 

{•£'t}i=l,...,258 = (J { E i}iel m ) ( 21 ) 
m=l,...,23 

with corresponding detector livetimes At m . The goal of this section is to include such 
time information, together with the set of monthly thermal reactor powers {-Pj^}, into 
a maximum likelihood analysis. The generalization is straightforward: monthly neutrino 
fluxes 4> m , signal spectra S m , background spectra B m , 6 and probability distributions D m are 
defined as 

d(f> m \ ^ \ - Pjm Of dNf 

dE u ~ f^ x f^ x \vL)E s dE v ' 1 ' 

S m (E)=enMAt m JdE u ^ J dE' M f E \ E>) r(E,E>) , (23) 

B m (E) = B(E) ^ , (24) 

S m (E) + B m (E) = N thco D m (E) , (25) 
respectively, fulfilling the relations 

J2S m (E) = S(E) , (26) 

rn 

Y.B m (E) = B(E) , (27) 

m 

and the probability normalization condition 

W dED m (E) = l. (28) 

m ^thr 



We prefer to plot the — currently small — allowed regions in linear scale, rather than in logarithmic scale. 
In particular, the log-scale in tan 2 6*12, introduced in 01 an d used in 0, can be usefully replaced by a 
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linear scale in sin 612 , which preserves the #i2_octant symmetry |4lj when applicable (this is not the case 

for a linear scale in tan 2 612, as used, e.g., in 18] ). 

We assume that all background components are constant in time. 
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The likelihood of the spectral shape information acquires then an explicit (monthly) time 
dependence, 

£sha P c = n n dm) , (29) 

m ielm 

while the functional forms of £ ra te and £ sys t remain the same as in Eqs. (fT5j) and (jTSjl . 
respectively. We have thus all the ingredients to calculate a likelihood function in energy 
(event-by-event) and time (monthly-binned). 

Notice that the likelihood function in energy and time reduces to the energy-only like- 
lihood function in the limit of constant reactor powers (-P^ = Pj h ), up to an irrelevant 
overall factor (the product of At m /At ratios); this limit provides a useful cross-check of the 
numerical results. 



A. Constraints on the oscillation parameters 

We start the discussion of the time- dependent effects in a case where they are (currently) 
not expected to play a significant role, namely, in the determination of the oscillation pa- 
rameters (5m 2 , sin 2 ^12). The KamLAND bounds on these two parameters are basically 
dominated by two different pieces of information: the energy spectrum shape and its nor- 
malization. In particular, the 5m 2 parameter governs the oscillation phase, which is strongly 
constrained by the observation of half-period of oscillations (7|,|33j|. This observation is still 
dominated by statistical errors 42], which currently hide subleading time-dependent effects, 
such as a possible shift of the "oscillation dip" for strong reactor power variations (as shown 
in Fig. 1). On the other hand, the sin 2 9 12 parameter governs the oscillation amplitude, 
whose bounds are dominated by normalization systematics 42 1. which are not reduced by 



adding time information. Therefore, within current uncertainties, we do not expect the 
mass-mixing bounds from the energy spectrum analysis (Fig. 3) to be significantly changed 
by adding time information. 

Figure 4 shows the results of our likelihood analysis in energy and time, which confirms 
the above expectations. A comparison with Fig. 3 reveals appreciable changes only in the 
"high-5m 2 " allowed region (so-called LMA-II solution [23|), which appears to be slightly 
more disfavored by adding time information. This trend allows to exclude with more confi- 
dence the LMA-II solution in combination with solar data (which, by themselves, still allow 
relatively high values of 5m 2 (18]). For the sake of completeness, and for later purposes, 
we show in Fi g. 5 the oscillation parameter bounds from our analysis of all current solar 
neutrino data |43j (including the latest full SNO spectral results (l8|) plus the KamLAND 
likelihood analysis in time and energy. The bounds in Fig. 5 contain, to our knowledge, 
the largest amount of solar and reactor neutrino information which is publicly available at 
present. 



B. Probing time variations of the reactor neutrino flux 

Variations of the reactor powers and of the livetime efficiency generate time variations 
of the event rate in KamLAND. Therefore, theoretical event rates including (not including) 
time information are expected to track more (less) faithfully the observed event rates. In 
Fig. 6 we plot the observed monthly counts in KamLAND, with respect to our calculated 



counts, 7 with and without reactor power variations. The comparison of the two panels shows 
at a glance that the correlation among the 23 points is more evident when monthly reactor 
powers (Pjm°) are included, with respect to the hypothetical case of constant reactor powers 
(Pjm = Pj h )- Quantitatively, the correlation index decreases from 0.73 (left panel) to 0.58 
(right panel). In the right panel, the correlation would be further reduced for hypothetically 
constant detector livetimes (At m = At/23), since all points would then collapse onto a single 
vertical line (not shown). 

The significant covariance between observed and calculated counts — when time informa- 
tion is fully included — suggests that KamLAND is indeed tracking reactor neutrino flux 
variations. In this sense, Figure 6 qualitatively agrees with the correlation test shown in the 
first figure of 0. We refrain, however, from fitting a "straight line" through the points in the 
left panel of Fig. 6, since we know of no clear way to include the point-by-point systematics 
and the large statistical fluctuations in such a linear fit. A maximum-likelihood test of time 
variations appears to be more appropriate, both to deal with small monthly counts and to 
include event-by-event energies and systematics. 

In order to test the null hypothesis of no time variations against the hypothesis of actual 
time variations of the reactor neutrino flux, we introduce an auxiliary variable r], interpolat- 
ing between the two cases. In this way, the hypothesis test is transformed into a parameter 
estimation test Formally, we assume that parameter rj modulates all reactor powers 

through the equation 

P*( V ) = Pf + r, AP^ , (30) 

where APj^ = Pjj^— P- are the actual power variations in each month. Thus rj continuously 
"switches on" reactor neutrino flux variations from the null case rj — (no time variations) 
to the real case rj — 1 (actual time variations). 

By using reactors powers P^Jjj) defined as in the above equation, we build a likelihood 
function in energy and time C(Sm 2 , sin 2 9u, rj), and marginalize it with respect to the oscil- 
lation parameters. The results are shown in Fig. 7, in terms of the function Ay^irf). The 
hypothetical case of constant averaged reactor powers (rj = 0) is definitely disfavored by 
KamLAND data, as compared with any case including time variations (0 < rj < 1). In par- 
ticular, the difference with respect to the case of actual time variations (77 = 1) amounts to 
about 3a (A% 2 ~ 9). We conclude that the results in Fig. 7 (and, to some extent, in Fig. 6) 
can be taken as a statistically significant indication that reactor neutrino flux variations 
have been seen in KamLAND. 

Finally, it is interesting to note that, in Fig. 7, the addition of solar neutrino information 
(through an additional A% 2 function which depends on (5m 2 , sin 2 6*12) but not on rj) does 
not significantly change the overall bounds on rj. In other words, as also observed in the 
previous subsection, energy and time information are largely decoupled in KamLAND (at 
present). The energy information indicates nonzero oscillation parameters, while the time 
information indicates nonzero variations of the reactor signal rate, with no appreciable cross- 
talk between these two pieces of information. Only with much smaller errors one might hope 
to see mixed effects (e.g., time-dependent changes of the energy spectrum dip). However, 
as we shall see in the next section, such "decoupling" of the oscillation parameters is not 
necessarily preserved in nonstandard cases, e.g., in a scenario with a hypothetical georeactor. 



7 Theoretical estimates refer to the solar+KamLAND best- fit oscillation parameters in Fig. 5. 



9 



IV. CONSTRAINTS ON THE GEOREACTOR 



It has been proposed [21| that there could be enough Uranium in the Earth's core to 
naturally start a nuclear fission chain over geological timescales, with a typical power (at 
the current epoch) of 3-10 TW jlJ], and possibly up to ~ 30 TW [45j]. The latter value is 
probably too high to be credible, since the addition of a typical radiogenic contribution of 
~ 20 TW (not to count other sources ji?}) would exceed the total Earth heat flux (esti- 
mated to be ~ 44 TW in 0] and recently revised down to ~ 31 TW in [49]). A georeactor 
power of ~ 10 TW is, however, comparable to the global Earth heat flux uncertainty |47j . 
and thus cannot be currently excluded by energy-budget arguments. On the other hand, 
there are independent geochemical and geophysical arguments which seem to disfavor any 
significant Uranium content in the core Hoi . Despite being largely ignored in the Earth 



science literature, the georeactor hypothesis has attracted some attention in the particle 
physics literature jHI, |5j| • 

In the KamLAND data analysis, a hypothetical georeactor can induce several effects. 
First, it increases the overall expected event rate. Second, it distorts the spectrum shape, 
both because its natural fuel composition can be significantly different from that of man- 
made reactors, and because the oscillation phase for L = i? e is different than for L ~ 0(100) 
km. Third, the georeactor signal is constant, while man-made reactors induce, in general, a 
variable signal in KamLAND. Therefore, we expect that a maximum likelihood analysis of 
the KamLAND data in energy and time, including the bounds on the oscillation parameters 
from solar neutrino data, can provide interesting constraints of the georeactor hypothesis. 
Technically, we implement the georeactor hypothesis by adding (in the KamLAND data 
analysis) a 21-th reactor at L = 6400 km, 8 with arbitrary constant power P geo . For definite- 
ness, we assume a current georeactor fuel ratio 235 U : 238 U ~ 0.75 : 0.25, with no significant 
Pu contribution |45j . 

Figure 8 shows the bounds on the oscillation parameters from our KamLAND maximum- 
likelihood analysis in energy and time, for the illustrative case P geo = 15 TW. The "wavy" 
contours of the lowest-5m 2 allowed region in Fig. 8 reflect the "ripples" created by georeactor 
neutrino oscillations on top the KamLAND energy spectrum (not shown). For the two 
allowed regions at higher values of 5m? , such (higher- frequency) ripples are smeared away by 
the finite KamLAND energy resolution, and the contours are smooth. More importantly, all 
the three allowed regions in Fig. 8 appear to be shifted to larger values of sin 2 9 12 , as compared 
with the standard (no georeactor) case in Fig. 4. This behavior is qualitatively expected, 
since larger mixing is needed to suppress the excess event rate due to the georeactor. 9 For 
increasing P geo , we should then expect an increasing tension with solar neutrino data, which 
fix sin #i2 around the value ~ 0.3 (as shown in Fig. 5) independently of P ge o- 

Let us now consider the results of a maximum-likelihood analysis where P geo is free, and 
the oscillation parameters are marginalized away. The results are shown in Fig. 9, in terms of 
the function Ax 2 (P geo ). From right to left, the four curves refer to increasingly informative 
and powerful analyses: (1) KamLAND likelihood in energy; (2) KamLAND likelihood in en- 
ergy and time; (3) KamLAND likelihood in energy, plus solar neutrino data; (4) KamLAND 
likelihood in energy and time, plus solar neutrino data. One can see that solar neutrino 
data can provide powerful (although indirect) constraints, by forbidding the large values of 

8 The georeactor radius is AL < O(10) km Q and thus negligible in this context (AL/L < 10~ 3 <C A_EZ£D. 

9 As a rule of thumb, a georeactor having power P geo = y TW increases the KamLAND rate by ~ y% |2£|. 
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sin 2 9 12 preferred for P geo > 0. To a lesser extent, the time information in KamLAND (consis- 
tent with known reactor source variations) also disfavor any additional constant georeactor 
contribution. In all four cases, we find no statistically significant evidence for P geo > 0, 
and can thus place meaningful upper bounds on its value. In particular, the most complete 
and powerful analysis in Fig. 9 (leftmost curve) provides the bound Pn- Pr , ^5 13 TW at 2a 
(95% C.L.), not too far from the typical expected range of a few TW j2l|. Basically, such 
bound tells us that, at ~ 2a level, the georeactor contribution should not exceed twice the 
KamLAND normalization uncertainty (i.e., ~ 13%). 

As a final remark we add that, since known reactor power variations help in constrain- 
ing a constant (hypothetical) georeactor neutrino flux, they can also be expected to help 
in constraining the constant (guaranteed |46|) geoneutrino flux below the current analysis 
threshold. In other words, as emphasized in [53j, a maximum likelihood analysis in both 
energy and time should provide a powerful tool for the statistical separation of the expected 
geoneutrino signal in KamLAND. Similarly, one might try to extend the current bounds on 
a (hypothetical) constant antineutrino flux from the Sun [H^ in the energy region where 
reactors provide a time-variable signal. 

Summarizing, we find that the inclusion of the (monthly-binned) time information in 
the KamLAND analysis corroborates the usual interpretation of the data, in terms of an 
oscillation-suppressed neutrino flux generated from known (time-variable) reactor sources. 
We find no indication for additional constant contribution from a natural georeactor, and 
place an upper limit P geo < 13 TW at 95% C.L. In any case, as emphasized in the Intro- 
duction, more refined and official KamLAND likelihood analyses (including, e.g., daily data 
about the detector and the reactors) will be crucial to improve and check such conclusions. 

V. CONCLUSIONS 

So far, published KamLAND data analyses have been focussed to the energy spectrum 
of neutrino events. In this work, after checking that we can reproduce in detail the official 
KamLAND likelihood analysis in energy, we have tried to add the time information to 
the analysis. In particular, by including monthly-binned data on Japanese reactor powers, 
KamLAND event-by-event energies, and detector livetimes, we find that the case of actual 
time variations of reactor powers is significantly preferred (~ 3a) over the hypothetical case 
of no time variations. This interesting indication is basically unaltered by adding solar 
neutrino constraints on the oscillation parameters. We have also considered the effect of a 
hypothetical georeactor with power P geo in the analysis. We find increasingly tighter upper 
bounds as more data (from time variations and from solar neutrinos) are included, down to 
P geo < 13 TW at 95% C.L. 

Our analysis supports the standard interpretation of the observed KamLAND neutrino 
events as generated by known reactors sources and affected by flavor oscillations with mass- 
mixing parameters consistent with solar neutrino data. Implications for a hypothetical 
constant neutrino source in the Earth interior (the georeactor) start to emerge. Other 
constant-flux sources (e.g., geoneutrinos or solar antineutrinos) might be usefully constrained 
in a similar way. We hope that these encouraging results may motivate other independent 
analyses of the reactor information in the time domain, especially by the KamLAND col- 
laboration that, by using the complete and fully controlled data set, can certainly provide 
more reliable results and explore interesting new facets of the topics touched in this work. 
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KamLAND positron energy spectrum 




FIG. 1: KamLAND absolute spectrum (without backgrounds) as a function of prompt energy. 
Solid curve: standard spectrum for best-fit LMA parameters. Dashed curve: spectrum with no 
contributions from the Kashiwazaki reactor power plant. Dotted curve: spectrum with additional 
contribution from a 15 TW georeactor. The vertical line indicates the analysis threshold (2.6 MeV). 
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Fission number flux 
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FIG. 2: Absolute fission number flux at KamLAND, as a function of the reactor distance. Dotted 
histogram: KamLAND estimate. Solid histogram: this work. 
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Likelihood analysis (energy) 
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FIG. 3: Bounds on the oscillation parameters from a maximum-likelihood analysis of the Kam- 
LAND energy spectrum (event-by-event). 
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Likelihood analysis (energy and time) 
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FIG. 4: Bounds on the oscillation parameters from a maximum-likelihood analysis of the Kam- 
LAND information in energy (event-by-event) and time (monthly binned). 
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Combined analysis (Solar + KamLAND) 
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FIG. 5: Bounds on the oscillation parameters from all current solar neutrino data and their 
combination with the KamLAND bounds in Fig. 4. 
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Correlation of observed vs calculated events in each month 
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FIG. 6: Monthly counts of events observed in KamLAND, plotted against the corresponding 
theoretical counts (calculated for the global best-fit LMA parameters in Fig. 5). Left panel: time 
variations of reactor powers included. Right panel: variations excluded (average powers used). 
The positive correlation between observed and calculated counts appears to be more pronounced 
in the first case. In both cases, monthly KamLAND livetimes are included. 
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FIG. 7: Bounds on the parameter 77, which interpolates between the extreme cases of no time 
variations of reactor powers (rj = 0) and actual time variations of reactor powers (rj = 1). The 
value rj = 1 is significantly preferred over r/ = 0. Bounds at 1, 2, and 3a can be obtained at 
A^ 2 = 1, 4, and 9 (dotted horizontal lines). The results are dominated by the maximum likelihood 
analysis of KamLAND data in energy and time (dashed line) and are not appreciably affected by 
including solar neutrino data. The oscillation parameters (5m , sin 2 9 12 ) are marginalized in both 
cases. 
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Likelihood analysis (energy+time) 
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FIG. 8: As in Fig. 4, but adding the contribution from a hypothetical georeactor with P geo = 15 
TW. 
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Bounds on georeactor power 
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FIG. 9: Bounds on the georeactor power from increasingly detailed analyses (with marginalized 
oscillation parameters). From right to left, the Ax 2 curves refer to: KamLAND analysis in energy; 
KamLAND analysis in energy and time; KamLAND analysis in energy plus solar neutrino data; 
KamLAND analysis in energy and time, plus solar neutrino data. At 95% C.L. (2a), the strongest 
upper bound (leftmost curve) is P geo ^$13 TW. 
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